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Abstract. Aims. We study the neutral hydrogen distribution and kinematics in the Local Group dwarf irregular galaxy IC 1613 
and compare them with the ionized gas distribution and stellar content of the galaxy. We discuss several mechanisms which 
may be responsible for the origin of the observed complicated HI morphology and compare parameters of the most prominent 
kpc-scale HI structure with the multiple SNe scenario. 

Methods. The observations were performed with the Vary Large Array of NRAO with a linear resolution ~ 23 pc at the adopted 
distance of 725 kpc and the spectral channel width of 2.57 km s"'. The numerical calculations have been provided with our 
2.5D Lagrangian scheme based on the thin layer approximation. 

Results. We found that the ISM of the galaxy is highly inhomogeneous and identified a number of intermediate-scale (200 pc - 
300 pc in diameter) HI arcs and shells having expansion velocities of 10 to 20 km s"'. Besides these shells, several giant holes 
and arc-shaped structures have been revealed, whose radii exceed several hundred parsecs. We found that parameters of the 
most prominent (Mhi = 2.8 x 10' Mq) kpc-scale structure and the level of the detected star formation activity are inconsistent 
with the multiple SNe hypothesis. 

Key words. Galaxies; irregular - Galaxies; individual - IC1613; ISM - kinematics and dynamics; ISM - supershells 



1. Introduction 

The neutral component of the ISM in star-forming dwarf and 
spiral galaxies presents a collection of holes, arcs, and shell- 
like structures (Davies & Tovmassian, 1963; Heiles, 1984; 
Brinks & Bajaja, 1986; Puche et al. 1992; Ehlerovaet al. 2004). 
Recent studies of late-type dwarf galaxies indicate that ~ 1/3 
of the galaxies present a high column density and sometimes 
broken rings with Nhi > 10^' cm^^ (g^jj ^ isr^^i 2002). The 
origin of HI holes and shell-like structures has been the sub- 
ject of debate for more than two decades, but still remains a 
controversial issue (see, for example. Block & Walter, 2005). 
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Several different mechanisms have been proposed for the 
formation of the observed HI distribution and kinematics in dif- 
ferent galaxies (see, for example, reviews by Tenorio-Tagle & 
Bodenheimer, 1988; van der Hulst, 1996). 

In the standard approach based on the Dyson & de Vries 
(1972), Weaver et al. (1977) model, HI shells result from the 
cumulative effects of multiple stellar winds and supernovae ex- 
plosions inside the volume presently encompassed by the shell. 
In such a case the energy released by stellar winds and super- 
novae is thermalized and generates a hot, high pressure cav- 
ity that drives a shock wave into the surrounding interstellar 
medium (ISM). This shock collects interstellar gas and com- 
presses it into a dense, expanding shell. As the hot cavity 
continues to expand, the temperature and pressure inside the 
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cavity drop while the outer shell decelerates. When the pres- 
sure becomes comparable to that of the ambient medium, the 
shock wave vanishes. At such a pressure-confined stage (Koo 
& McKee, 1992), the shell stalls and starts to disintegrate with 
its local sound speed. 

Despite many cases in which this approach gives a plausi- 
ble explanation for the observed structures (see, for example, 
Oey & Garcia-Segura, 2004; Vorobyov et al., 2004, 2005; van 
der Hulst, 1996 and references therein), it was recognized long 
ago (Heiles, 1984; Tenorio-Tagle & Bodenheimer, 1988) that 
this scenario cannot explain all types of observed structures. In 
particular, it encounters problems when attempting to explain 
the origin of the largest supershells where the stellar associa- 
tions found within the HI rings seem to be inadequate energy 
sources. More recently, Rhode et al. (1999) found that in the 
case of the Holmberg II (HoII) galaxy, the observed upper lim- 
its for stellar cluster remnants detected inside neutral hydrogen 
holes is in many cases inconsistent with the SNe hypothesis. 
Although a year later Stewart et al. (2000) re-examined their 
results using deep far-ultraviolet (FUV), Hor and HI data and 
decided that the energy deposited by stellar winds and SNe ex- 
plosions into the ISM of HoII is sufficient, there are many HI 
holes that have no optical counterparts consistent with the rate 
of mechanical energy required by the standard model. 

In other galaxies, Kim et al. (1999) found only a weak 
correlation between the LMC neutral hydrogen holes and the 
HII regions and concluded that the hypothesis of multiple 
winds and supernovae is inconsistent with their data. Perna 
& Gaensler (2004) proposed a non-traditional approach to the 
problem by counting the number of radio pulsars associated 
with Milky Way supershells. They performed Monte-Carlo 
simulations for the expected pulsar population and predicted 
that several tens of radio pulsars may be detected in the direc- 
tions to the largest Milky Way supershells, and few of them 
should be found inside regions encompassed by the contours 
of smaller shells. They found that for the smaller shells their 
predictions are consistent with the detected pulsar population. 
However for the giant supershell GSH242-03-I-37 the observed 
number of radio pulsars is inconsistent with the multiple winds 
and supernovae model. 

Several other mechanisms that do not require violent stellar 
activity include collisions of high velocity clouds with galactic 
disks (Tenorio-Tagle, 1981), the radiation pressure from field 
stars (Elmegreen & Chiang, 1982), the non-linear evolution of 
self-gravitating turbulent galactic disks (Wada et al. 2000; Dib 
& Burkert, 2004), the ram pressure of the intergalactic medium 
(Bureau & Carignan, 2002), and even exotic mechanisms such 
as the distortion of the ISM by powerful gamma-ray bursts 
(Efremov et al. 1999). Despite many efforts attempting to re- 
solve this problem, the detailed comparison of different sce- 
narios with objects detected in different galaxies remain scarce 
and controversial. 

IC 1613 is a member of the Local Group and is located 725 
- 730 kpc from the Milky Way (Freedman, 1988a,b; Dolphin 
et al., 2001). Because of its proximity , it represents an excel- 
lent target to study the interplay among the stellar and gaseous 
components in galaxies. Its HI mass, derived from single-dish 
observations by Lake & SkiUman (1989), is 6 X 10^ Mq. The 



galactic rotation is slow, with an almost linear rotation curve 
inside the inner 2.5 kpc radius, having a maximum amplitude 
of rotation velocity around 25 km s"' (Lake & Skillman, 1989). 

The stellar population of IC 1613 and its star formation his- 
tory (SFH) have been studied by a number of authors. Hunter 
et al. (1993) estimated the total Ha luminosity to be Lho ~ 
3 X 10^** erg s The associated star formation rate (SFR) is in 
good agreement with that obtained by Cole et al. (1999), who 
found SFR~ 3 x 10"^ Mq yr ' over the past 10 Myr. Skillman 
et al. (2003) did not find any evidence for recent episodes of vi- 
olent star formation in this galaxy from their analysis of color- 
magnitude diagrams and concluded that it looks unlikely that 
this and similar galaxies may be responsible for the chemical 
enrichment of the intergalactic medium. Georgiev et al. (1999) 
obtained deep UBV photometry for more than 3000 stars in se- 
lected areas of IC 1613. The isochrone fit of their data clearly 
indicates the presence of young (age < 30 Myr) stars and stellar 
associations. 

On the other hand, Fabry-Perot interferometric studies 
of the ionized gas kinematics (Valdez-Gutierrez et al. 2001; 
Lozinskaya et al. 2003) revealed a number of supershells cov- 
ering the whole extent of the optical galaxy. The complex ra- 
dial velocity pattern indicates expansion velocities up to 60 - 
75 km s The majority of the supershells presents an inte- 
rior OB-association, suggesting a physical link between them. 
Lozinskaya et al. (2003) identified three large (~ 300 pc), 
bright HI supershells in the north-east part of the galaxy. They 
expand at ~ 15 km s ' and overlap with several smaller HII 
shells distributed along the walls of the HI supershells and are 
ionized by the young OB-associations. The location of the HII 
shells and OB-associations suggests that sequential star forma- 
tion was triggered in this region by larger neutral hydrogen su- 
pershells. This implies that, despite a quiet history of star for- 
mation, the ISM of IC 1613 has complex gas kinematics and 
morphology which are Unked to the embedded stellar popula- 
tion. 

Here we present a study of the neutral hydrogen distribu- 
tion and kinematics within this galaxy based on Very Large 
Array observations and compare them with our previous stud- 
ies of the ionized gas and stellar content in this galaxy. We 
have found several giant HI rings and arc-like structures. The 
smallest ones present clear 21 -cm line splitting indicating that 
in these cases we observe expanding shells generated by the 
embedded OB-associations. The origin of the largest ones is 
not clear. In particular, the largest, ~ 1 kpc in diameter, struc- 
ture contains several OB-associations of different ages which 
appear to be potential sources of energy. However we did not 
find any indication of a regular shell expansion that leads to an 
assumption that this structure may represent a case of the stand- 
ing, pressure-confined shell. This initiated us to compare the 
observed parameters of the "main supershell" with the standard 
model based on the multiple winds and supernovae hypothesis. 
We found that the SNe hypothesis is in poor agreement with the 
observationaUy restricted rate of star formation that indicates a 
different origin of this large, kpc-scale structure. 

The paper is organized as follows. In Section 2 we present 
the results of the observations, the procedures for the data re- 
duction, discuss the neutral hydrogen distribution and the ob- 
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served velocity field and compare this with properties of the 
ionized ISM and stellar content of the galaxy. In particular, 
we discuss properties of the largest HI structure 1 kpc in di- 
ameter that is similar to the largest HI supershells detected in 
other galaxies. In Section 3 we discuss the standard superbub- 
ble model, present the results of numerical calculations and 
compare them with the observed parameters of the "main su- 
persheU", the largest HI structure detected in IC 1613. Section 
4 summarizes our results. 

2. HI data: morphology, kinematics and tlieir linlcs 
with other galactic components 

Neutral hydrogen observations were made with the Very Large 
Array (VLA) of The National Radio Astronomy Observatory 
(NRAO) in configurations B, C and D. The application for the 
current observations was compiled by E.M. Wilcots. The first 
results, together with the brief discussion for data reduction, 
were presented in Lozinskaya et al. (2001). 

The data cube has an angular resolution of 7" .4x7" .0, and a 
pixel size of 3" x 3". This implies a linear resolution of ~ 23 pc 
at the adopted distance of 725 kpc. The spectral channel width 
is 2.57 km s~'. Our data cube allows us to derive the distri- 
bution and kinematics of the HI gas of the entire galaxy. To 
build up the position-velocity diagrams and the integrated 21- 
cm line intensity map, the 40 velocity channels, from the 127 
available, have been selected by the condition that the signal- 
to-noise ratio exceeds the 2cr limit in three or more adjacent 
velocity channels (a windowing technique; see, for example, 
Thornley & Mundy, 1997). The selected channels cover veloc- 
ity interval from -279 to -178 km s"^. 

2. 1. The distribution of neutral liydrogen 

The integrated 21 -cm line flux density derived from our data 
is 482Jykms"'. This implies that the total HI mass of the 
galaxy is 6.0 x 10^ M©. This value is in good agreement with 
the single-dish measurements by Lake & Skillman (1989; see 
Table II from their paper), despite a "negative bowl" in our 
data that leads to a negligible underestimation of the HI mass. 
An average ID velocity dispersion of 5.9 km s ' was mea- 
sured from our data cube after accounting for the instrumen- 
tal resolution. This corresponds to a 3D velocity dispersion 
CIS M — 10.2 km s"' . Lake & Skillman (1989) adopted a some- 
what larger value of 7.5 km s"' for the ID velocity dispersion. 
However, our data have better spectral and spatial resolutions 
(5 km s"' versus 6.2 km s"', and 7" versus 45"). 

The map of the 21 -cm Une surface brightness is presented 
in Figure la. The three brightest and most clearly defined shells 
(Shells I, II and III) in the N-E sector of the galaxy have diam- 
eters of 200 - 300 pc and are clearly related to the only region 
of noticeably enhanced star formation discovered so far in the 
galaxy (see Lozinskaya et al. 2003). Several more shells and 
arcs of similar sizes having lower surface brightness are clearly 
identified in Figure la and Figure 3a (e.g. Shells IV and V). 

Besides these intermediate- size shells, one can easily rec- 
ognize in Figure 1 several larger arcs and ring-like structures 
which are similar in size and morphology to the supergiant 



Table 1. The distribution of the HI mass in the galaxy. 
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shells revealed by Meabum (1980) at the Large Magellanic 
Cloud. The most prominent structure seen in Figure la co- 
incides with the giant region of low surface brightness emis- 
sion visible in the lower angular resolution map of Lake and 
Skillman (1989). It is a hole surrounded by a thick HI ridge. 
The hole has a 5' diameter (~ 1 kpc diameter at the distance 
adopted for IC 1613) and is centered at RA(2000) = Ih 4m 52s, 
Dec(2000) = 2°07'. The thickness of the HI ridge, AR, is about 
200 pc. The Northern, Eastern and Southern parts of the ring 
have a complicated morphology and a kinematical pattern sug- 
gesting a feedback from the local star formation. Shells II and 
III form the brightest part of the northern wall. The northwest 
sector of the ring appears to have a lower surface brightness 
and is not clearly identified in our map. Figure la also presents 
the locations of the OB-associations (shown by circles whose 
radii are proportional to the sizes of the OB-associations) taken 
from the list of Borissova et al. (2004). Figure lb displays a 
population of superbubbles identified by Valdez-Gutierrez et 
al. (2001). Hereafter the whole structure is referred to as the 
"main supershell". 

We spUt the galaxy into a number of segments (Figure Ic) 
and then calculated the mass of each segment by integration of 
the HI column density over the segment surface. The segment 
masses are presented in Table 1. We associate the "main super- 
shell" with segments 1, 2, 11, 14, 17 (Figure Ic). The HI mass 
of the "main supershell" is then 2.8 x 10^ M©, which corre- 
sponds to a total mass of 3.9 x 10^ Mg when including helium. 

Another example of a large shell-like structure is an ellipti- 
cal hole in the HI distribution at the north of the main supershell 
surrounded by sectors 1, 5, 6, 7, 9 and 10 (Figure Ic). One can 
also recognize several incomplete shells formed by sectors 10, 
1 1 and 12, sectors 14, 1 1, 12, 13 and 15, and sectors 13, 15 and 
16 at the northwest and west of the main supershell. The size of 
the shell to the north is similar to that of the "main supershell" 
but it has a smaller mass. Mm ~ 9.7 x 10'' Mq. 

2.2. The ionized component, the stellar population, 
and the local SFR 

The examination of the distribution of the HII regions and su- 
perbubbles from the list of Valdez-Gutierrez et al. (2001; Tables 
4 and 5) shows that many of them are located inside the area 
encompassed by the contour of the "main supershell" or di- 
rectly inside the HI ridge. Figure lb displays the boundaries of 
the superbubbles found by Valdez-Gutierrez et al. (2001) su- 
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Table 2. OB-associations from Borissova et al. (2004) that are 
located within the field of the main supershell 



10 20 30 40 50 60 70 80 90 100 

M (MJ 



Fig. 2. The local star formation rate. The number of massive 
stars, members of OB-associations, which are located inside 
the main supershell (solid Une) as compared to the number 
of massive stars derived from the constant star formation rate 
model (dotted line). The reasonable fit to the high mass top of 
the observed mass function requires the local rate of star for- 
mation to be SFRiocai = 4.5 x 10"'' Mq yr"i. 

perimposed on the HI image of the galaxy. Several superbub- 
bles (Rll, R12, R14, R15, R16 and R17) are located inside the 
"main supershell". Rl, R4 and R6 seem to coincide with HI 
shells I, n and III identified by Lozinskaya et al. (2003). The 
rest are projected onto the HI ridge and seem to be associated 
with the extended region of star formation at the northeast of 
the galaxy. 

The Ha luminosities of the Hll regions and superbubbles 
obtained by Valdez-Gutierrez et al. (2001) allow us to estimate 
the mass of the ionized gas and the local star formation rate 
within the area encircled by the "main supershell". The mass 
of the ionized gas associated with superbubbles was calculated 
from the relation Mhii = 3.3 x 10"^ 2;(/?,/lpc)^ (M©), 
where radii of superbubbles, and electron number densities, 

were taken from Table 8 of Valdez-Gutierrez et al. (2001), 
the thickness of the superbubble shells was assumed to be 1/12 
of their radii, and the contribution from the helium was taken 
into account. For the HII regions we assume a homogeneous 
density distribution. Then Mmi = 3.3 X 10"' Yj(Ri/^Pcf ^ej. 
The rms electron number densities were calculated from the 
relation n^j = (2.3 X lO'^F/fa,,/©^/?,)!/^^ where is the equiv- 
alent angular diameter of the Hll region expressed in square 
arc-seconds and Fnaj is the Ha flux in the units of ergs cm"^ 
s"'. The radii of the superbubbles and the Hll regions were 
measured in parsecs. The total mass of the ionized component, 
Mhii = 1.6 X 10* M©, is small in comparison with the HI mass 
of the main supershell and, therefore, we neglected its contri- 
bution in our calculations. 

We then estimated the local SFR using the caUbration given 
inKennicutt(1998) 



SFRiM^yr-') = 7.9 X lO~'*'LHa(erg s'') 



42, 
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adding up the luminosities of the superbubbles (Rll, R12, R14, 
R15, R16 and R17 from the list of Valdez-Gutierrez et al. 



(2001), see Figure lb) that are located inside the inner bound- 
ary of the "main supershell" and therefore may contribute to 
the dynamics of the supershell. We found that the redden- 
ing correction is small because, according to the NASA/IPAC- 
Extragalactic-Database (NED), is only 0.083 mag. The Ha 
luminosity then is 3.24 x 10^' erg s and the local star forma- 
tion rate, SFRiocai - 2.6 x 10 Mq yr This value has been 
compared with the SFR (Figure 2) obtained from the analysis 
of stellar contents of OB-associations from the list of Borissova 
et al. (2004). The masses, the absolute visual magnitudes, Mv, 
and the ages of OB-associations within the field of the main 
supershell are Usted in Table 2. 

The ages of the associations have been calculated from the 
observed color-magnitude diagrams which were converted into 
absolute magnitudes and the true colors for each of the associ- 
ations. The foreground reddening was calculated for each star 
individually using the reddening free Q parameter and the lat- 
est calibration equations of Massey et al. (2000) summarized 
in their Table 4. A distance modulus and a metallicity were 
adopted where (m - M)o = 24.27 ± 0.1 (Dolphin et al. 2000), 
and Z = 0.004 (Mateo 1998). The ages were then calculated by 
the best fit to the theoretical isochrones from the Geneva library 
(see Lejeune & Schaerer 2001). The Ukely error of our method 
is 2-2.5 Myr. 
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In order to obtain the masses of individual stars and a gen- 
eralized mass function we have used the procedure proposed by 
Massey et al. (2000), which allows the determination of the ef- 
fective temperatures (Tgif) and the bolometric corrections (BC) 
for individual stars. The evolutionary tracks from Lejeune & 
Schaerer (2001) for Z = 0.004 were then used to establish mass 
intervals. The local SFR derived from the comparison of this 
mass function with the constant star formation rate model is 
SFRiacai^ 4.5x10-^ Mo yr-\ 

Obviously, these values only give a lower limit to the star 
formation rate within the galactic area covered at the current 
epoch by the "main supershell". The mass function obtained 
from the Borissova et al. (2004) data is incomplete at the low 
mass end and does not account for the field stars. The local SFR 
derived from the Ha data is also related to the current epoch 
and does not consider the diffuse Ha emission. Nevertheless 
our estimates are reasonable when compared to the mean star 
formation rate over the galaxy. Therefore we adopt them as a 
likely limit to the SFR within the field of the main supershell. 

2.3. The HI kinematics 

In order to study the kinematics of the interstellar medium 
(ISM) in IC 1613 we obtained position-velocity diagrams in 
more than 30 directions. The orientations of several scans and 
the corresponding position-velocity diagrams are presented in 
Figure 3. 

The kinematics of Shells I, II and III were discussed by 
Lozinskaya et al. (2003). Shells II and III have expansion ve- 
locities of 12 - 18 km s"' and kinematical ages 5.3 - 5.6 Myrs. 
The expansion velocity of Shell I could not be measured, while 
its size is similar to that of Shells II and III. Shells IV and V 
mentioned in the previous section are much fainter and their 
properties are harder to determine than those of Shells I, II and 
III. Their radii are about 150 - 200 pc. The position-velocity 
diagram taken across the eastern part of the "main supershell" 
(cut 4 at Figure 3a) reveals an offset of approximately of 15 
km s~' between the regular rotation velocity of the galactic gas 
(which is in between -225 km s"' and -235 km s"' in this re- 
gion) and the maximum values found in the field of the shells. 
This is -215 km s"' for shell IV (at 170" - 230" along the cut), 
and -220 km s"' at the field of Shell V (at 250" - 315" along 
the cut). These offsets together with the characteristic shapes of 
local velocity ellipses suggest that the expansion velocities of 
both of the shells are around 15 km s ' . 

The position-velocity diagram taken along a shell formed 
by sectors 10, 1 1 and 12 also indicates a regular shell expansion 
(see cuts 7 and 8 at Figure 3, at 550" - 700" and 300" - 400" 
along the cuts). The expansion velocity is around 13 km s 
The size of this shell is 3' that corresponds to a linear diameter, 
D = 640 pc, and the HI mass is 6.4 x lO*" Mq (see Table 1). 
Other giant structures at the north and northwest of the "main 
supershell" show no evidences of regular shell expansion. 

The "main supershell" was cut in 12 directions. We did 
not find any indication of regular expansion of the "main su- 
pershell" in any of the cuts. Only local gas perturbations were 
identified in some sectors of the "main supershell" (e.g. at the 



positions corresponding to 115" - 200" of cut 2, 135" - 180" 
of cut 3, 430" - 630" of cut 5 and 90" - 210" of cut 7). The 
observed local velocities are in the range 215 - 245 km s"^ 
The best example of the local high velocity flow is observed 
in the vicinity of the WO star which is located at the southeast 
wall of the "main supershell" (see Lozinskaya et al., 2001). In 
this region we probably observe the blow-out of the WO stellar 
wind from the shell into the low density bubble interior. The 
local velocities observed in this region (see Figure 3b, at the 
positions corresponding to 190" - 300" along the cut 1) fall 
into the range -205 km s"' - -250 km s"' . We believe that other 
arc-shaped structures seen in the inner part of Segment 17 re- 
sult from the collective feedback given by the associations 26, 
34, 36, 37, 38, 42, 45 of Borissova et al. (2004), which are lo- 
cated near the ridge of the "main supershell". The expansion 
velocities of the local arcs are less than or equal to 25 km s"' . 

The absence of regular expansion and the large thickness 
of the "main supershell" led us to conclude that if this struc- 
ture results from the combined action of the embedded groups 
of massive stars, then the main supershell has already reached 
its standing position and is now confined by the external gas 
pressure (see Koo & McKee, 1992). 

3. The multiple SNe scenario 

In the standard SNe scenario, the shell expansion is supported 
by a high thermal pressure provided by multiple supernova ex- 
plosions and stellar winds from massive stars recently formed 
inside the superbubble volume. In the case of an homogeneous 
ISM and a constant input rate of mechanical energy, the time 
evolution of the radius, R, the expansion velocity, V, and the 
mass, M, of the shell are given by (see Bisnovatyi-Kogan & 
Silich, 1995, and references therein): 



R 



V 



iU5\ 
\l54nl 
3R 
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47T 



1/5 



'-'mech 
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M = —R'pisM-- 



(2) 

(3) 
(4) 



where t is the bubble age, L,„ec/7 is the rate of the mechanical 
energy deposited by SNe explosions and stellar winds and pisM 
is the density of the ISM. Eliminating the evolutionary time t 
and the ISM density pisM from equations - (0), one can 
find the mechanical power, L„ech, required to create a shell with 
radius R, mass M and expansion velocity V: 
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100 



10 -1 

erg s 



(5) 



where M-j is the mass of the shell in units of 10^ Mq, Viq is the 
expansion velocity in units of 10 km s"' and Z^ioo is the radius 
of the shell in units of 100 pc. 

The transition to the pressure-confined stage occurs when 
the velocity of the shock approaches the ambient gas turbulent 
speed, V - ctjsm- At that time the shell begins to expand in 
both directions from the stand-off radius, Rstand, with the lo- 
cal sound speed, c.,/,. This may be consistent with the irregular 
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boundary of the main supershell and with the standard scenario. 
If this is the case, it takes the shell ~ 0.5AR/Csh > lO' yr to 
reach the observed thickness. Thus the detected configuration 
may be observed for a good fraction (~ 1/4) of the total evo- 
lutionary time, t. The stand-ofF radius is Rmner < Rstmd < ^out^ 
where Rimer is the irmer and Rgm is the outer radius of the shell. 
We assume that 

A/? 

R.tand = Rinmr + — = 600 /?C , (6) 

where A/? ^ 200 pc is the mean thickness of the shell. For the 
rate of mechanical energy we then obtain: Lmech « 1.9 x 10^^ 
erg s"'. This energy input rate requires a coeval star cluster 
of approximately lO'^Mg with a Salpeter initial mass func- 
tion with a lower and upper mass Umits of 1 Mq and IOOMq. 
We do not see stellar remnants of such massive clusters in the 
field of view of the "main supershell". On the contrary, OB- 
associations found inside the ring have a wide range of ages 
(see Table 2) that indicates that in the region encompassed by 
the main supershell the star formation lasts for more than 30 
Myr. This suggests that a continuous star formation rate would 
be a better approximation for the star formation activity that 
occurred in this region. 

In this case the rate of mechanical energy cannot be ap- 
proximated by a constant value during the whole evolution. In 
such a scenario the maximum value is reached ~ 30 Myr af- 
ter the beginning of star formation (see Leitherer & Heckman, 
1995). On the other hand the thermal pressure of the interstellar 
gas also reduces the shell radius below the standard value (Oey 
& Garcia-Segura, 2004). The calculated energy may be also af- 
fected by the deviation of the ISM density from a homogeneous 
distribution. 

To take into consideration these efl'ects and to obtain a re- 
alistic estimate for a SFR that is required to build up a shell 
of similar to the observed structure, we have provided sev- 
eral numerical runs using for calculations our 2.5D Lagrangian 
scheme based on the thin layer approximation (see Bisnovatyi- 
Kogan & Silich 1995; Sihch & Tenorio-Tagle 1998). As an in- 
put model we adopt an exponential galactic disk with a constant 
velocity dispersion a = 10.2 km s"\ 

piSM = Po exp {-\z\IH,) , (7) 

where po is the midplane ISM density. It was assumed that the 
disk is embedded in a low density gaseous halo whose density 
was taken to be tihaio - 10"^ cm"-'. We neglected the small shell 
distortion (Palous et al., 1990; Silich, 1992) provided by the 
slow galactic rotation. The Z-component of the gravitational 
field was calculated from the equilibrium conditions. The mass 
and energy deposition rates were calculated from our simph- 
fied starburst model (see Silich et al. 2002) assuming that the 
SFR remains constant until the shell waist attains the observed 
radius. The SFR and the midplane ISM density po were then 
calculated by iterations from the condition that the expansion 
velocity of the shell waist becomes equal to the ISM turbulent 
velocity and the calculated mass of the shell approaches the ob- 
served value when the radius of the shell waist reaches 600 pc. 

The calculations have been provided for two different (thin 
and thick) galactic disk models with the characteristic scale- 
heights of the exponential disk = 200 pc and = 500 



pc, respectively. Figure 4 presents the evolutionary sequence 
of superbubble shapes for these two cases. 

In both cases the bubbles grow faster in the direction of the 
steepest density gradient. In the case of the thin disk, the su- 
perbubble already reaches a typical hour-glass form after 20 
Myr (Figure 4a). Soon the shell accelerates rapidly into the 
low-density galactic halo and the top of the shell is disrupted 
under the action of Rayleigh-Taylor instabiUties. At this age 
(« 20 Myr), the top of the shell reaches approximately three 
characteristic Z-heights of the interstellar gas distribution. The 
remnant is then able to "blow out", driving its high tempera- 
ture gas between the fragments into the halo of the host galaxy 
(Silich & Tenorio-Tagle, 2001), and to establish the cone-like 
structure seen in Figure 4b. The bubble then loses its thermal 
pressure and the waist of the shell reaches its standing position 
at the age t ^ 34 Myr. Despite the large size of the cone, the 
bulk of the swept-up material is located within a thin cylindri- 
cal segment around the base of the cone. In this particular case, 
about 50% of the mass of the remnant is concentrated inside 
a 400 pc (200 pc above and 200 pc below the midplane of the 
galaxy) segment. 

To fit the observed parameters of the main supershell the 
midplane gas number density should be no « 2.8 cm"^ and the 
star formation rate SFRw 7.5 x lO^-'M© yr"'. This value is 
more then an order of magnitude larger than the local star for- 
mation rate, 5 F/Jfoca/ « (3-4)xl0"'' Moyr"', derived from the 
analysis of our Ha data and stellar population of the embedded 
OB-associations, and approximately three times larger than the 
average star formation rate obtained by Cole et al. (1999) for 
the entire galaxy. Thus the thin-disk model is completely in- 
consistent with the the standard SNe scenario. 

In the thick disk case the shell remains almost round after 
20 Myr of expansion. Even at the end of the calculations, at 
t ^ 37Myr, the shell is only slightly distorted by the ISM den- 
sity gradient (see Figure 4d). The Rayleigh-Taylor instabilities 
never disrupt the top of the shell and the hot gas remains bound 
within the cylindrically-shaped remnant slowly expanding in 
the Z-direction. To fit the observed parameters (the stand-off 
radius and the mass of the shell), the mid-plane density of the 
disk should be * 1 .4 cm"-' and the star formation rate must be 
about 2.3 x 10"^ M© yr "^ This value significantly (~ 5 times) 
exceeds our limit to the local star formation rate and is only 
slightly below the mean star formation rate obtained for the 
entire galaxy. There are also no indications on a massive sin- 
gle stellar cluster that may support shell expansion down to the 
observed radius. This leads us to the conclusion that multiple 
SNe hypothesis is inconsistent with the observed parameters of 
the "main supershell" also in the case of the thick disk ISM 
distribution. Some different mechanism (e.g. the outward force 
initiated by the radiation pressure of the field stars, Elmegreen 
& Chiang, 1982) is required to build up the observed structure. 

4. Conclusions 

Here we presented the results of our high resolution study of 
neutral hydrogen distribution and kinematics in a Local Group 
member, the dwarf irregular galaxy IC 1613 and compared 
them to our previous examination of the stellar population and 
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ionized component of the ISM in this galaxy. The analysis of 
the 21 cm line data cube reveals a complicated, highly inho- 
mogeneous structure of the interstellar medium with a number 
of HI holes, shells and arc-like structures of different sizes and 
expansion velocities. 

We found that smaller shells with diameters of a few hun- 
dred parsecs tend to expand with velocities lOkms"' < y < 20 
km s"' whereas the largest structures do not present evidences 
of global expansion. The only exception is the large incomplete 
shell with mass ~ 6.4x 10'' Mq and expansion velocity, V ~ 13 
km s ' at the northwest part of the galaxy. The best example of 
the large standing structure is an HI hole 1 kpc in diameter sur- 
rounded by a 2.8 x 10^ M© neutral hydrogen ring centered at 
RA(2000)=lh4m 52s, D(2000)=2deg 07'. 

There are many OB-associations of different ages inside the 
hole whose energetics potentially may be responsible for the 
development of this structure. To examine this hypothesis we 
have compared the energy required by the multiple supernovae 
model with that inferred from the star formation rate derived 
from the analysis of the embedded stellar population and H„ 
emission. We found that parameters of the observed structure 
and the level of the star formation activity detected inside the 
region encompassed by the contour of the "main supershell" 
are inconsistent with the SNe hypothesis. This indicates that, 
as it has been revealed in the past for the largest HI holes in 
HoII and LMC (Rhode et al. 1999; Broun et al. 1997), some 
different mechanism is required to build up the largest kpc-size 
HI structure that is detected in IC 1613. 
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Fig. 1. The distribution of neutral hydrogen. Upper left: the 21-cm Une brightness map with OB associations from the list of 
Borissova et al. (2004) superposed. OB associations are shown by circles whose radii are proportional to their sizes. The coor- 
dinates and radii of associations are taken from the Table 1 of Borissova et al. (2004). Upper right: the 21-cm hne brightness 
map of the southeast sector of the galaxy with the superbubbles of Valdez-Gutierrez et al. (2001) superposed. Lower panel: the 
distribution of the HI column density throughout the galaxy split into a number of segments. In all panels north is at the top and 
east is on the left. 
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Fig. 3. The neutral gas kinematics. Upper panel displays orientations of cuts used to study the neutral gas kinematics. Cuts are 
identified by numbers. The region of active star formation in the northeast sector of the galaxy is shown in the upper left-hand 
box. The position of the WO star is marked by a cross. Panels at the bottom represent the position- velocity diagrams along each 
cut indicated at the upper panel. All velocities are heliocentric. 
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Fig. 4. The shape of the supershell. Panels a), b) and c), d) present the results of the calculations for the thin {H^ = 200 pc) and 
thick (i/j. — 500 pc) disk models, respectively. The ages of the shells in panels a) and c) are similar 20Myr). Panels b) and d) 
present supershells at the moment when they reach their standing positions at the midplane of the galaxy. 



